Nuclear Magnetic Relaxation Rate in a Noncentrosymmetric Superconductor 
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For a noncentrosymmetric superconductor such as CePtsSi, we consider a Cooper pairing model 
with a two-component order parameter composed of spin-singlet and spin-triplet pairing compo- 
nents. We demonstrate that such a model on a qualitative level accounts for experimentally observed 
features of the temperature dependence of the nuclear spin-lattice relaxation rate , namely a 
peak just below T c and a line-node gap behavior at low temperatures. 



Inversion symmetry is one of the key points for the 
formation of Cooper pairs in superconductors. Unusual 
properties arise in superconductors whose crystal struc- 
ture does not possess an inversion center (e.g., Refs. 
jH HI El, d Q i an d references therein) . The recent dis- 
covery of superconductivity in the noncentrosymmetric 
heavy Fermion compound CePtsSi has initiated much in- 
terest, as experimental data revealed various intriguing 
features The essential element in modelling non- 

centrosymmetric systems is the presence of antisymmet- 
ric spin-orbit coupling |3, [10( . One of the characteristic 
features in the superconducting phase is the mixing of 
the spin-singlet and spin-triplet Cooper pairing channels 
which arc otherwise distinguished by parity [2(|. This is 
likely responsible for the surprisingly high value of the 
upper critical field H C 2 which exceeds the paramagnetic 
limit HE1ISIIJE3- CePt 3 Si displays further intriguing 
properties. Recent nuclear magnetic resonance (NMR) 
experiments found an overall anomalous temperature de- 
endence of the nuclear spin-lattice relaxation rate T^ 1 
" EEl The behavior of Tf 1 shows a (Hebel-Slichter) 
peak just below the superconducting critical temperature 
T c , and simultaneously a T 3 dependence at low tempera- 
tures indicating line nodes in the quasiparticle gap. Such 
a gap with line nodes is also suggested by measurements 
of the London penetration depth At first sight, those 
experimental results seem to be mutually contradicting, 
as the features of unconventional Cooper pairing (line 
nodes) and of conventional superconductivity (peak in 
T^ 1 due to the coherence effect) are implied at the same 
time by the temperature dependence of T-, -1 @, E10]- 

In the present study, we demonstrate that these appar- 
ently conflicting behaviors can be reconciled by taking 
account of the mixing of the pairing channels with op- 
posite parity, which naturally occurs in superconductors 
without inversion center. For this purpose, we consider 
a pairing model with a two-component order parameter 
consisting of spin-singlet and spin-triplet pairing compo- 
nents. Such a model contains the necessary ingredients 
to account for the observed features of Tj -1 , i.e., the line- 
node behavior at low temperatures and the coherence ef- 
fect just below T c . At the same time, the pairing model 



would also explain the temperature dependence of the 
London penetration depth qualitatively and be consis- 
tent with earlier studies of H C 2 j3> E| • 

We base our analysis on a Hamiltonian considered in 
Ref . , in which the lack of inversion symmetry is incor- 
orated through the antisymmetric spin-orbit coupling 
12, 16] 



agk ■ S with 



9k = \/^{-kyi k oi 
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Here, S is the electron spin operator, a (> 0) denotes the 
strength of the spin-orbit coupling, the vector (g-k = 
—gk) is determined by symmetry arguments, and k = 
(k x ,k v ,k z ) = (cos <f> sin 9, sin <f> sin 9, cos 9). 

Starting from the tetragonal symmetry of CePtaSi, it is 
possible to classify the basic pairing states, distinguishing 
the spin-singlet and spin-triplet states 0, . The gen- 
eral argument by Anderson 17] shows that the inversion 
symmetry is a key element for the realization of spin- 
triplet pairing. Hence, turning on antisymmetric spin- 
orbit coupling would be detrimental for spin-triplet pair- 
ing. The detailed examination of this aspect, however, 
led to the conclusion that among the spin-triplet states, 
only the d vector parallel to gt, [ dk \\ (—k y , k x ,0) ] , is 
robust against this symmetry reduction effect of the spin- 
orbit coupling |3|]. Moreover, the antisymmetric spin- 
orbit coupling mixes spin-singlet and spin-triplet pair- 
ings. Interestingly, the conventional s-wave spin-singlet 
pairing state mixes precisely with the state corresponding 
to dk || gk }16j |. which will be essential for our discussion. 

We consider a superconducting gap function in the fol- 
lowing mixing form, 

A fe = (*(T + d k ■ &)ia y 



*<T 



(2) 



with the s-wave pairing component \& and the d vector 
dk = A(—k y , k x , 0). While the spin-triplet part has point 
nodes, this pairing state can possess line nodes in a gap as 
a result of the combination with the s-wave component. 
In this paper, we adopt the isotropic s-wave pairing as \P 
for simplicity. 

We will discuss the superconducting phase by means of 
the quasiclassical theory of superconductivity p^.Hfjll2rj| . 
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The quasiclassical Green function g is written as a matrix 
in Nambu (particle-hole) space, 



g(r, k,iuj n ) = -in 



9 ^ if 
-if -9 j 



(3) 



The vector r is the real-space coordinates, the unit vec- 
tor fe indicates the position on the Fermi surface, and 
u> n = irT(2n + 1) is the Matsubara frequency. Through- 
out the paper, "hat" • denotes the 2x2 matrix in the 
spin space, and "check" • denotes the 4x4 matrix com- 
posed of the 2x2 Nambu space and the 2x2 spin space. 
The Eilenberger equation which includes the spin-orbit 
coupling [Eq. Q] is given by [U H El H 

iv F ■ V<? + [mj„t 3 - ag k ■ S - A fe , g] =0, (4) 
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Here, <x = (a Xl <7 y ,<T z ) is the Pauli matrix, vp is the Fermi 
velocity, and [a,b] — ab — ba. We use units in which 
h = /cb = E The Eilenberger equation is s upp lemented 
by the normalization condition, g 2 = —ir 2 l flftlElT]. 

For the pairing state [Eq. (J2J1 , we obtain the following 
Green functions in the case of spatially uniform system, 

9 = Sioi + ffii^ii, (9a) 

/ = (Mi + fn&iijivy, 



f 



/(Ml + /Mil) 



9 = —&y(9lfa+Sii<rii)<ry> 
with the matrices 01,11 defined by 0, 0, 

9k = ( 
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Here, k = (k X} h y} 0) — (cos 0, sin 0,0), and 
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(9c) 
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(10) 
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* + Asin0 f-Asintf 
/i = „ > /n = 5 > ( llb ) 



Si 



$* + A*sin0 - **-A*sin0 , . 

h = 5 > /n = 5 . ( llc ) 



9l 



9u 



The denominators -Blii are given as 



Si,n = v^ + l*±Asin0| 2 , 



(lid) 



(12) 



and the signs in front of the square root are determined 
by the condition, sgn(Re{giji}) = sgn(Re{o;„}). Note 
that the Green functions [Eq. Qll[)1 do not explicitly de- 
pend on the spin-orbit coupling constant a. This result 
of the Eilenberger equation reflects the fact that the spin- 
triplet component contained in the pairing state [Eq. J2J] 
is not affected by the antisymmetric spin-orbit coupling 
of Eq. © 0). 

The above Green functions labeled by the indices I, 
II belong to the two distinct Fermi surfaces which are 
split by the spin-de gen eracy lifting due to the spin-orbit 
coupling agk ■ S 0,13 Hij]- While in general the density 
of states on the two Fermi surfaces is different, we assume 
that the difference is small and ignore it here. 

The pairing interaction leading to the gap function 
[Eq. 0] 

can be characterized by three coupling con- 
stants, A s , At, and A m . Here, A s and At result from the 
pairing interaction within each spin channel (s: singlet, 
t: triplet). X m appears as a scattering of Cooper pairs 
between the two channels, which is allowed in a system 
without inversion symmetry |l6| . The gap equations are 
written as 



X s ttT (f+) + Kn*T ]T (sin0/_), (13) 



A = A t 7rT J2 (sm0/-) + A m7 rT £ (/+), (14) 



\uj n I <u; c 



\uj n | <UJ C 



where f± = (/i ± fu) /2, and u c is the cutoff energy. The 
brackets (• ■ ■ ) denote the average over the Fermi surface. 
We assume here the spherical Fermi surface for simplicity. 
In the limit T — > T c , the linearized gap equations allow 
us to determine At and X s by 



At = 



v\ T , 

2\w 



w 



\T r . 



2 

A. s = - At + 



E 

0<n<(w o /wT-l)/2 



AV J 



2n + 1 



(15) 



(16) 



(17) 



Bi 



Bu 



if the parameters X m and v are given. 

In Fig. 1, we show the temperature dependence of the 
order parameters A and \& obtained from the gap equa- 
tions [Eqs. 03J) and ifTIj)]. We set u c = 20T C , X m = 0.12, 
and v = 0.5, yielding A t w 0.39 and A s ~ 0.16. When A 
is fixed to be real and positive without loss of generality, 
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FIG. 1: The temperature dependence of the order parameters 
in units of T c . The spin-triplet component A (solid line) and 
the spin-singlet s-wave component ^ (dashed line). lu c — 
20T C , A m = 0.12, and v = 0.5. Both A and * are real and 
positive. 



the solution such that 'F is also real and positive is stable 
for the above parameters. 

We consider now the nuclear spin-lattice relaxation 
rate divided by the temperature, 1/T±T. Assuming for 
the NMR experiment a static magnetic field in the z di- 
rection, we arrive at the following expression for \jT\T 
in terms of the quasiclassical Green function g [2(| , 



Ti(T c )T c 
T X (T)T 



1 

4T 



coslr (w/2T) 



W(u), (18) 



with 



W{w) = <a 22 »>(a^(-u;)) (af T M><a^(- W )),(19) 
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g(r, fe, iuj n — > uj + irf) 
- g(r,k,iuj n -> u - irj) , (20) 
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FIG. 2: The temperature dependence of the nuclear spin- 
lattice relaxation rate 1/T\T (solid lines). r\ — 10 _4 T C . (a) 
Dashed line is the contribution of the anomalous Green func- 
tions Wff related to the coherence effect. Dash-dotted line is 
the contribution of the regular Green functions Wgg related 
to the density of states, (b) Plot of the same data on a double- 
logarithmic scale. Dotted line is a plot of T 2 . From the plot, 
it is noticed that T^ 1 follows the T 3 law at low temperatures. 



where r\ (> 0) is an infmitesimally small constant, and 
we set t] = 1(T 4 T C . 

The temperature dependence of 1/T\T is shown in Fig. 
2(a). To obtain it, we used the temperature dependence 
of A and \& shown in Fig. 1. Obviously, \fT\T (solid 
line) possesses a peak just below T c . In order to iden- 
tify the origin of this peak, in Fig. 2(a) we also plot the 
contributions of the two terms in Eq. (jT§|) separately: 
W = W GG + W FF , 



W GG (u) = (af l (oj))(a\\(-u ) )), 
Wff(u) = -{a 2 l \(iu)){a] 2 l (-u;)). 



(22) 



(23) 



W GG and Wff are composed of the regular Green func- 
tions and the anomalous Green functions, respectively. 
The coherence factor is represented as 1 + Wff /W gg . 
The contribution of Wff (dashed line) is related to the 
coherence effect and gives the dominant contribution to 
the peak below T c . In contrast, W GG (dash-dotted line) 
describes the contribution of the density of states. The 
contribution to the peak from the singularity of the den- 
sity of states at the gap edge is minor, since this singu- 
larity is rather weak due to the anisotropy of the gaps 
| ^ ± A sin 8 1 on both Fermi surfaces I, II. Note that the 
anisotropy in the real compound would be enhanced be- 
cause of the anisotropy of the Fermi surfaces . There- 
fore, the peak in the total \jT\T (solid line) can clearly be 
attributed to the coherence-factor-induced enhancement 
of the relaxation rate T^ 1 originated from the coherence 
effect. 

Turning to the low-temperature behavior, we present 
the same data on a double-logarithmic scale in Fig. 2(b). 
The temperature dependence of 1/TjT exhibits a T 2 
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FIG. 3: Schematic figures of the gap structure on the Fermi 
surface, (a) Cross section of the Fermi surface, showing the 
places of gap nodes at which the gap |*t/ — Asin^l = 0. (b) 
Line nodes of the gap in the three dimensional k space. 

power law at low temperatures, characteristic of the pres- 
ence of line nodes in the gap. These nodes are the result 
from the superposition of spin-singlet and spin-triplet 
contributions (each separately would not produce line 
nodes). On the Fermi surface I, the gap is | M/ -f- Asin^l 
and is nodeless, referring to Eq. Ullbjl . (note that *Jy > 0, 
A > 0, and < 9 < tt). On the other hand, it has the 
form — A sin 9\ on the Fermi surface II, where nodes 
can appear (see Fig. 3). For *Jy| <C |A| the gap nodes 
are located at small 9 close to the poles on the Fermi 
surface, leading to nearly point nodes. With growing 
s-wave component *]/, however, the gap nodes move to- 
wards the equator of the Fermi surface and form line 
nodes if 1*51 < |A|, and they disappear if |*F| > |A|. 
These line nodes on the Fermi surface II lead to the low- 
temperature T 3 law in Tf 1 (i.e., T 2 in 1/TiT) as shown 
in Fig. 2(b), which is in qualitative agreement with the 
experimental result. Moreover, the power law observed 
in the London penetration depth [Sj is also consistent 
with this explanation. 

In conclusion, we identify a natural Cooper pairing 
state [Eq. (J2J] consisting of a spin-singlet s-wave *P and 
a spin-triplet A components. This pairing state explains 
the set of presently available experiments consistently. 
They include the early experiments on the upper critical 
field H C 2 whose behavior can be well explained on the 
basis of this pairing state 0, E, El- The line nodes 
which can occur due to the superposition of the two 
spin channels are very compatible with the observation 
of the power law behaviors at low temperatures in the 
London penetration depth and the NMR relaxation rate 
Ti 1 @, US E3- In addition, this state allows for the 
Hebel-Slichter coherence peak in Tj -1 . It is unlikely that 
the observed peak could originate from a singularity of 
the density of states at the gap edge, since such fea- 
tures are most likely washed out by the anisotropy of 
the quasiparticle gap in the real material. It should also 
be noted that in CePtsSi the superconductivity coexists 
with an antiferromagnetic phase. Muon spin relaxation 
experiments show that there is little mutual influence of 
such two orders, suggesting that those two orders might 



be associated with different Fermi surfaces respectively 
d, ^| . This aspect has to be taken into account when 
low temperature thermodynamics is analyzed in this ma- 
terial. However, the London penetration depth contains 
exclusively the information on the superconductivity, and 
therefore measurements of it Q belong to the cleanest ex- 
periments in probing the gap topology. This fact gives 
confidence in the existence of line nodes in CePtsSi. 
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Note added. — After completing this work, we noticed 
that a similar idea concerning the coexistence between 
the coherence effect in T^ 1 and the line node in a gap was 
briefly discussed in Ref. [27j by Fujimoto very recently. 
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